Scattering of unstable particles in a finite volume: 
the case of up scattering and the ax (1260) resonance 
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We present a way to evaluate the scattering of unstable particles quantized in a finite volume with 
the aim of extracting physical observables for infinite volume from lattice data. We illustrate the 
method with the scattering which generates dynamically the axial-vector ai(1260) resonance. 
Energy levels in a finite box are evaluated both considering the p as a stable and unstable resonance 
and we find significant differences between both cases. We discuss how to solve the problem to 
get the physical scattering amplitudes in the infinite volume, and hence phase shifts, from possible 
lattice results on energy levels quantized inside a finite box. 



I. INTRODUCTION 

The determination of hadron spectra is one of the chal- 
lenging tasks of Lattice QCD and much effort is being 
devoted to this problem [ll4l2||. Large pion masses are 
commonly used in these calculations [j, [l3l - H7| . The 
"avoided level crossing" is usually taken as a signal of 
a resonance, but this criteria has been shown insufficient 
for resonances with a large width [l8l - [20| . A more accu- 
rate method consists on the use of Liischer's approach, 
for resonances with one decay channel, in order to pro- 
duce phase shifts for the decay channel from the discrete 
energy levels in the box [2l|, [22j • This method has been 
recently improved [2(| by keeping the full relativistic two 
body propagator (Liischer's approach keeps the imagi- 
nary part of this propagator exactly but makes approx- 
imations on the real part) and extending the method to 
two or more coupled channels. The new method also 
combines conceptual and technical simplicity and serves 
as a guideline for future lattice calculations. Follow ups 
of this new practical method have been done in (23| for 
the application of the Jiilich approach to meson baryon 
interaction and in [24| for the interaction of the DK and 
r]D s system where the D s *o(2317) resonance is dynam- 
ical ly g enerated from the interaction of these particles 
[2a - l2q | . The case of the k resonance in the Kir channel 
is also addressed along the lines of [2(| in [29| . 

The case of scattering of unstable particles deserves 
a special care since in the box one must also discretize 
the momenta of the decay products of all the particles. 
One such system would be the 7rA system where the A 
is allowed to decay into ttN. The generalization of the 
work of [2(| to this problem has been done in [3(| • 

The problem of scattering of unstable particles will 
have to be faced by the lattice QCD calculations. So 
far, problems which would require this treatment have 
been studied assuming stable particles. This is the case 
of the irp scattering, from where the ai(1260) resonance 
is qualitatively obtained, assuming the p to be a stable 
particle, in an actual lattice QCD simulation using the 
first two levels for a fixed size of the box [31|. In the 



present paper we face directly this problem and provide 
the formalism to address it, also for the case of an unsta- 
ble p resonance. For this purpose recall that in the chiral 
unitary approach the axial vector ai(1260) resonance is 
dynamically generated from the interaction of irp and 
KK* in coupled channels [3i3. l33|. where np is the dom- 
inant channel. Then we follow the approach of Ref. [33| 
and solve the interaction of 7rp within a box with pe- 
riodic boundary conditions. These boundary conditions 
are imposed on the spectator it and on the two it that 
come from the p meson decay. For the it — 2tt system 
we choose the global center of mass (CM) frame, but the 
two pions that lead to the p in the 7T7t loop function are 
in a moving frame and this forces us to make the dis- 
cretization of the levels in this moving frame, a problem 
which is well studied in [34]; [H[ . Furthermore the p is 
a p-wave resonance which requires a different method to 
obtain the selfenergy than the s-wave resonances. All 
these problems will be dealt with in the present paper. 



II. FORMALISM 

In the chiral unitary approach the scattering matrix in 
coupled channels is given by the Bethe-Salpeter equation 
in its factorized form 



T= [l-VGl^V = [V~ 



(1) 



where V is the matrix for the transition potentials be- 
tween the channels and G is a diagonal matrix with the 
i th element, Gi, given by the loop function of two meson 
propagators, a pseudoscalar and a vector meson, which 
is defined as 



Gi 



1 



1 



d 4 p 

(2tt) 4 (P - pf - Mf + ie p 2 - mf + 



, (2) 



where m, and Mi are the masses of the two mesons and P 
the four-momentum of the global meson-meson system. 
Note that in Eq. (J2]) we have not considered the possi- 
ble widths of the mesons, thus it is only valid for stable 
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mesons. The unstable mesons case will be addressed in 
section IIVI 

For the case of scattering of a pseudoscalar with a vec- 
tor meson, for instance Tip, KK*, as in the present case, 
the interaction is taken from the chiral Lagrangians (36j 
and has the form 1331 



V PV 



V 



(3) 



where e, e ' are the polarization vectors of the initial and 
final vector mesons. The term e ■ e ' factorizes in all terms 
of the Bethe-Salpeter series, V, VGV, etc., and finally 
in the T matrix. Hence, we omit this factor in what fol- 
lows. The explicit expression of the potentials, properly 
projected onto s-wave, is thus [33j 



"8f y 
-\(Mf 



3s - (Mf + m\ + M] + m)) 



m?)(M, 2 -m|) 



(4) 



where / = 92.5 MeV is the pion decay constant, the index 
i(j) represents the initial (final) PV state in the isospin 
basis and Mi(Mj) and mi(rrij) correspond to the masses 
of the initial (final) vector mesons and initial (final) pseu- 
doscalar mesons, for which we use an average value for 
each isospin multiplet. The explicit values of the numeri- 
cal coefficients, Cy , can be found in Ref. [13]. For the np 

isovector amplitude, which we need for the present work, 
r i=i _ _ 9 

The loop function in Eq. ^ needs to be regularized 
and this can be accomplished either with dimensional 
regularization or with a three-momentum cutoff. The 
equivalence of both methods was shown in Refs. nam. 
In dimensional regularization the integral of Eq . (12"]) is 
evaluated and gives for meson-meson systems [371 l39j 



Gj(s, m,-, Mi) 



(4tt) 2 



cii(p) + loj 



Mf 



mf + s 



2s 



log 



Mf 



log (s - [Mf - m 2 ) + 2^Q,(^)) 



+ log (s + (Mf - mj) + 2V^Qi(Vs)) 

- log (-s + (Mf - mf) + 2yTsQi(Vs)) 

- log (s - (Mf - m\) + 2V5Qi(Vi)) 



(5) 



where s — E 2 , with E the energy of the system in the 
center of mass frame, Qi the on shell momentum of the 
particles in the channel i, (i a regularization scale and 
ai(p) a subtraction constant (note that there is only one 
degree of freedom, not two independent parameters). 



In other works one uses regularization with a cutoff in 
three momentum once the p° integration is analytically 
performed and one gets 



G, 



d 3 p 



1 



(p)+uj 2 (p) 



W\<Piy 



(2tt) 3 2wi(p) u 2 (p) £ 2 - ( Wl (p) + u 2 (p)Y + ie 



(6) 



with mi , m 2 corresponding to vtii and Mi of Eq. @ . 

When one wants to obtain the energy levels in the fi- 
nite box, instead of integrating over the energy states of 
the continuum with p being a continuous variable as in 
Eq. © , one must sum over the discrete momenta allowed 
in a finite box of side L with periodic boundary condi- 
tions. We then have to replace G by G = diag (G\, G 2 ), 
where 



, |p|<Pm i 
P 

2tt 



1 



UJl(p) + UJ 2 (p) 



2uji(p) u 2 (p) E 2 ~ (ui(p) + uj 2 (p)) 2 



p = — n, n£Z J 



(7) 



This is the procedure followed in [20(. The eigenener- 
gies of the box correspond to energies that produce poles 
in the T matrix, Eq. (fTJ), which correspond to zeros of 
the determinant of 1 — VG. 



III. ONE CHANNEL ANALYSIS 

In the present problem the threshold of KK* is above 
the mass of the ai(1260) resonance and, as found in [33l |. 
the 7rp channel is more important than the KK* one. 
For this reason we shall perform the analysis with just 
the irp channel, as also done for the lattice calculation of 

a 

The one channel problem can be easily solved and is 
very simple, as shown in (20| . The T matrix for infi- 
nite volume can be obtained for the energies which are 
eigenvalues of the box by 



T(E) = (V~HE)-G(E)Y 



G(E) - G(E) 



(8) 

since G(E) = V^ 1 (E) is the condition for the T matrix 
to have a pole for the finite box. 
Hence we find 



T(E)- 1 = lim 

p max — >oo 



P<Pa 



(9) 



where I(p) is the integrand of Eq. ^ 
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Hp) 



2ui(p) uj 2 (p) E 2 - (uii(p) +uj 2 (p)) 2 + ie 



(10) 



This result is the one obtained in Ref. [20| starting with 
cutoff rcgularization and, as proved in Ref. [201 ] . it is noth- 
ing else than Liischer formula [2l|,|22|, except that Eq. (|9|) 
keeps all the terms of the relativistic two body propaga- 
tor, while in Liischer's approach one neglects terms in 
Re I(p) which are exponentially suppressed in the physi- 
cal region, but can become sizable below threshold, or in 
other cases when small volumes are used or large energies 
are involved. 



IV. GENERALIZATION TO SCATTERING OF 
UNSTABLE PARTICLES 



P 2 = (P-P«) 2 



Pi 



with Mr 



2Pp„ =s + ml- 2 v ^E(p) 
(12) 



/s7 the invariant mass of the two pion sys- 



tem. We have chosen the CM for the irp system and 
hence p v = —p. 

Since II(p) is a Lorentz invariant magnitude one can 
evaluate it in the CM frame of the p meson. However, 
the analogous magnitude II (p) in the finite box must take 
into account the boundary conditions for the it momenta 
in the moving frame. 

In order to calculate the p-mr vertex, let us consider 
the standard Lagrangian for the coupling of one vector 
to two pseudoscalars |4l[ 



In this section we extend the approach to the case 
where we have one unstable particle. We shall work with 
the case of one channel, but the generalization to many 
coupled channels is straightforward. The consideration 
of the unstable particle requires to reevaluate the loop 
function G of Eq. ^ by using the dressed meson propa- 
gator including its selfenergy that accounts for the decay 
channels. This means substituting 



1 



u( P y 



(11) 



where H(p) is the meson selfenergy of the unstable parti- 
cle. Diagrammatically this means that we must evaluate 
the loop diagram of Fig. [TJ 



p\ P 



q 



p-q 



FIG. 1. The Tip loop diagram considering the p meson selfen- 
ergy. 



In order to calculate this loop function we must first 
evaluate the p selfenergy, II (p). One knows [32|,|33[ that 
differences between the sum and the integral in Eq. §§§ in 
regions of p far away from the pole of I(p) are exponen- 
tially suppressed in L. The sizeable differences stem from 
regions of p close to the on shell momentum p on where 
I(Pon) has a pole. For this reason we evaluate U(p) for 
values of p where the tt and irir systems can be placed 
on shell. This is the same prescription taken in Ref. [30] | . 
Then we have 



C 



VPP 



-i 5 y([P,0 M P]V> 



(13) 



where gy = M p /2f 7r , with f n = 92.5 MeV the pion de- 
cay constant, P, V*, the ST/ (3) matrices of the pseu- 
doscalar and vector mesons and (...) standing for the 
SU(3) trace. From this Lagrangian we find 



t, 



2V2 



gvq-e, 



(14) 



with q taken in the p meson rest frame and e the p po- 
larization vector. The p selfenergy is then given by 



*n(p) = 



d 4 q 



(2tt) 4 q 2 - ml + ie (P - q) 2 - ml + ie 
)2V2gvq-e(-i)2V2g v q-e'. (15) 



x(- 

In the p rest frame, where we evaluate the pirn vertex, we 
have qieiqje'j, and for symmetry reasons we can replace 
qiqj by q ' £y/3. Furthermore, in the chiral unitary ap- 
proach it is justified to use the on-shell approach where 
the vertices are factorized by their on shell form. Thus 
we finally obtain 



(16) 



where |<f on | = y/ s p /A — ml and Gt 77T (s p ) is the loop func- 
tion of Eq. (|2|) for two pions. The function G V7T (s p ) can 
be regularized by means of a cutoff, <7 max , in the mod- 
ulus of the three-momentum q and its explicit analytic 
expression is [38] 



nog 



8 2 -1 1 



a In ■ 



or 



1 



or — 1 



2 In 



"(I 



(17) 



2 

max' 



where a = \/ l — Amj/s n and r = y/T+rn%/q- 

In Ref. [42| the function is evaluated in dimen- 
sional regularization. This is equivalent to removing the 
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divergent part of II (s p ) of Eq. (fT?)) and substituting it 
by a subtraction constant which is then constrained by 
experimental data. In Ref. [421] , where the p meson is 
studied within the chiral unitary approach, one finds that 



CT+ 1 



G^(s p ) = abx-^-+b 



a - 1 



(18) 



with 

b = -2 + d\ = -2- 



'■k 



, mi 1 , m K 1 

ln / + 2 ln 



p 2 2 
(19) 

and p — m p . One way to get the result of Eq. (fIS|) is to 
take the limit 



G»„(s p )= lim 



(4tt) 2 



2 In 



2<z„ 



(20) 

since 2 ln(2g max /m 7r ) cancels the divergent part of the 
square bracket in Eq. (fT?]) when q max — > oo. 

When we evaluate the selfenergy in the finite box, 
IT(p), we shall use this expression but G^x will be re- 
placed by the corresponding discrete sum, G^tt- For 
the numerical evaluation, the quantity inside the square 
bracket in the previous equation has to be evaluated for 
high values of g max in order to get the convergence. How- 
ever, it oscillates around the convergence value for not 
very large values of qWx- Hence an average over the 
interval q max ~ [2,2.8] GeV numerically gets the conver- 
gence value and thus we take this average in the numer- 
ical evaluation, (see the analogous and further explained 
reasoning in Ref. (24l|). 



DISCRETIZATION OF THE p SELFENERGY 
IN THE MOVING FRAME 



The p selfenergy in the finite box, H(p), is given by 
Eq. (|16[) . where now G W7r will be given by Eq. (|2"0"|) evalu- 
ating G 7r7r (s p ) as a function of <? m ax for the finite box. In 
order to perform this evaluation we must substitute the 
d 3 q integration implicit in the evaluation of G Tr7r (s p ) in 
the infinite volume by a discrete sum over the pion mo- 
menta allowed in the finite box with appropriate bound- 
ary conditions which we take to be periodic. The subtlety 
is that the momenta q on in Eq. (|16[) and the variable of 
integration in Eq. (|15[) are in the CM of the two pions 
and the boundary conditions are in the box, where the 
pair of pions move with total momentum p. By perform- 
ing the q° integral in Eq. (|15|) we find in infinite space 
that 



n(p) = 8gU* 



^2 



d\ 



(2tt)= 



(21) 



l<?*l<9n 



with p 2 



and q* the momentum in the CM of the 



two pions. 

We must write the boost transformation from q to q*. 
By applying the Lorentz transformation from a moving 
frame with momentum p to a frame where the tttt system 
is at rest 1431 we find 



Q 1,2 = 9l,2 



Mr 



- 1 



Ql,2-P Ql,2 



Mr 



P- 



(22) 



where Mf = s p = p 2 , p° 2 = Mf + p 2 and the subindexes 
1,2, represent the two pions of the decay of the p meson. 
Demanding that q{ + q% =0 enforces + Q2 = P° or 
equivalently qf* + q^ = Mj. We take for q^°, q'^ , the on 
shell pion energy for the decay of an object of mass Mi 
at rest into two pions 



*o 

^1,2 



Mf 



'1,2 



'2,1 



2M T 



(23) 



Only this prescription makes ql° = q^ for two pions as 
we should expect for two identical particles in the center 
of mass. This provides then the boost for the off shell 
momenta in the loop, where q is arbitrary but the energy 
is the on shell one. Since we need the Jacobian of this 
transformation, it is useful to write Eq. (f2"2"|) in terms of 
the CM energy of the pion and we find 



1 1,2 = 91,2 



Mi 



91,2 -P 



*0 
9l 2 



p. (24) 



This equation is the one used in [301 ]. Furthermore we 

where the factor 



^jjf is the Jacobian of the transformation, with q the 
7r momentum in the tttt moving frame and then J Tg^yg 

becomes ^ - in the box. In summary, we must do the 
substitution 



must substitute J t^^U by J 



d^q Mj 



dh 



(2tt) 3 



1 ^M, 

8 



2tt. 



n, n £ 



(25) 



for the evaluation of the selfenergy in the box. 

When summing over q, the integrand takes the follow- 
ing structure for symmetry reasons: 



= eie'JaSij + bpiPj) 



(26) 



Eq. (|2"6")) is exact for any vector p placed along any of the 
axis and quite accurate, although not exact for vectors 
in other directions. Yet, the term of b is quite small, 
since the scale is 0((p/m p ) 2 ). Such terms have been 
systematically neglected in the approach of [33| and so 
do we here. By contracting Eq. (|26p after removing eie'j 
with Sij on one hand, and with PiPj on the other hand, 
we get two equations from where we find 



5 



« = ^/(r,g)^ 2 (i-cos 2 0), 

where cos 9 = p ■ q*/\p\\q*\- Hence we finally get 

l9*|<9m, 



(27) 



„ 1 - - < 



Mr 



1 — COS 



pVr(g*) s p -4ul(q*) 



up to the trivial e* • e' factor. 

Finally, in order to get a pole of the p propagator in 
the infinite volume at the physical p mass we subtract 
to n(sp) and n(s p ) the selfenergy Ren(m p ). Thus we 
replace 



n( Sp ) ^n( Sp )-Ren( m 2), 

U(s p ) -^U(s p )-BeU(m 2 X 



(29) 



VI. INCLUSION OF THE p SELFENERGY IN 
THE FINITE BOX AND INFINITE VOLUME 



We come back to the original problem of the irp inter- 
action with the p dressed by its selfenergy. In the approx- 
imation that we did to write I± p for the case where the 
three intermediate pions are placed on shell, the variable 
s p depends on p, see Eq. (TT2"]) . Hence, the p a integration 
of Eq. ^ when we replace the vector meson propagator 
as in Eq. pip can be performed in the same way as with 
two free propagators and we obtain 



G« P (E) 



d 3 p 



|p|<Pn 



(2tt) 3 2wio;2 



£2-( Wl + c, 2 )2-^n( Sp ) 



(30) 



where ui\ = uj\(p) and o>2 = ^(p) are the 7r and p on shell 
energies, and we have kept the leading term in n/wf, 
which is a small quantity close to the p on-shell from 
where the finite volume corrections stem mostly. 
For the box we substitute G^ p by G„- p given by 



_ ^ \P\<Pa- 

Gnp{E) = — 



1 



2wiW2 



^- (wi + W2 )2_^n( Sp ) 



(31) 



In one channel the scattering matrix in infinite volume 
is given by 



T 



v- 



G 



(32) 



and for the finite box 
T 



y- 



(33) 



The eigenenergies of the unstable irp system in the box 
are given by the energies that satisfy 



V- 1 = 



(34) 




(,oo 



1000 1200 

E [MeV] 



1600 



up 



FIG. 2. Loop function in the box G np (solid line) and V 1 
(dashed line) for L — 2m~ 1 and p max = 1 GeV. The dashed 
dotted line corresponds to the case with stable p, n(s p )=0. 



In Fig. [21 we plot V 1 and G Tp as a function of E for 

L = 2m~ 1 and p max = 1 GeV. The G np plot is very dif- 
ferent to a typical meson-meson loop function in infinite 
volume. It shows clear poles coming from zeros in the 
denominator of Eq. (|3"Tj) and from poles of n(s p ). These 
poles of G^p are not present in the infinite volume since 
in the integration the poles of the integrand provide an 
imaginary part to the loop function but not a pole af- 
ter performing the integration. However, this is not the 
case in the finite box since we do not have an integral 
but a summation. The intersection between both plots 
of Fig. [2]just provides the 7rp scattering eigenenergies in 
the box. 

It is interesting to note that the spectra obtained is 
also qualitatively different to the one obtained for the 
stable p which would be given by the intersection of V^ 1 
and G wp for the case of the stable p (dashed-dotted line 

in the figure). We can see that for stable p one has G,r P 
going to infinity when the energy E approaches one of the 
free energies of the irp system in the box. Then V^ 1 cuts 
G^p only once in between two neighboring asymptotes. 
When we discretize the tttt system, n(s p ) becomes infi- 
nite for the discrete energies of the moving tttt system in 
the box. With n(s p ) becoming infinite, with plus and mi- 
nus sign, in the denominator of Eq. (|31[) . independently 
of the value of E 2 — (oji +UJ2) 2 , close to the pole of n(s p ) 
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there will be an energy where the denominator will van- 
ish, leading to a pole of G^ p . Thus, we get asymptotes of 
G^p for values of E close to the free eigenenergies of up 
and also close to the free eigenenergies of the moving tttv. 
We observe that in between two asymptotes correspond- 
ing to the irp free eigenenergies (dashed dotted lines) one 
new asymptote has appeared corresponding to a free irir 
eigenenergy of the irir moving frame. As a result of it, 
the line V~ x cuts now two lines corresponding to the un- 
stable p and only one corresponding to the stable p. If we 
go to higher energies we observe that in between the next 
two asymptotes corresponding to the free eigenenergies 
of np there are now three extra asymptotes correspond- 
ing to free eigenenergies of 7T7T in the moving frame. This 
leads now to four eigenenergies of the interacting irp sys- 
tem in the box when we cut these lines with V . There 
is, thus, a proliferation of eigenenergies as a consequence 
of considering the p as an unstable state. 



VII. RESULTS 




800 - 



L [m 1 ] 

FIG. 3. The first six energy levels as a function of the cu- 
bic box size L for stable p-meson (II = 0) (dashed lines) and 
for unstable p-meson (solid lines) using p max = 1 GeV. The 
dotted lines indicate the free up energies of the box for com- 
parison. 



A. Energy levels in the box 

We first show in Fig. [3] the solutions of Eq. (|34|) which 
represent the energy levels for different values of the cubic 
size, L. The solid lines represent the full calculation, con- 
sidering the unstable p-meson via the consideration of its 
selfenergy in the box, II in Eq. (f3"Tj) . and the dashed lines 
are the same calculation but setting II = in Eq. (|3"Tj) . 
ie., considering the case of stable p meson. There is an 
infinite number of levels and we have plotted the first six 
ones for illustration. This distribution of energy levels 
depends on the cutoff p max and in particular this figure 
has been evaluated using p max = 1 GeV which is of the 
order of the cutoff used in Ref. [33j to get dynamically 
the axial- vector resonances. 

It is worth noting the significant difference between 
the consideration of the unstable p meson in comparison 
to the stable case. The main difference is that the levels 
tend to decrease and shrink as L increases and that there 
are extra levels between those already present in the sta- 
ble case, for instance between the second and third levels. 
This latter feature comes from the appearance of extra 
poles in G^ p due to the poles in II and this produces 
extra intersections with the smooth function V^ 1 as ex- 
plained above. This different behavior between the stable 
and unstable analysis must be considered as a caveat for 
lattice calculations that only consider the stable case. 



B. Inverse problem: getting np amplitudes and 
phase shifts from lattice-like data 

By "inverse problem" we refer to the problem of getting 
the actual scattering amplitudes (and hence by-product 
magnitudes like phase shifts) in the infinite space from 



data consisting of points over the energy levels in the 
box in the E vs. L plots, which is what a lattice calcu- 
lation would provide. In our case we can "synthetically" 
simulate this lattice-like data from our model generating 
points in the levels of Fig. [31 We will call this data lattice- 
like set although in the present work it is generated from 
our model for illustrative purposes. 

In Fig. 0] we represent by error bars the set generated 
for a particul ar ele ction of L values for the stable, 4(a)| 
and unstable, |4(b)| p cases, to which we have assigned a 
reasonable error of 10 MeV. The meaning of the shad- 
owed error bands will be explained later in this section 
when explaining the fit method pol ] , hence they must be 
forgotten for the moment. We are aware that it is difficult 
for a present lattice calculation to get such quantity of 
points as considered in Fig. [51 but the method is equally 
valid for less points and we have chosen such a set just 
to illustrate more clearly the method. In an actual in- 
verse problem the lattice-like generated set would just be 
replaced by actual lattice results. 

In Refs. (§3, HSHH Hil several methods were suggested 
to solve the inverse problem, mostly based on fitting the 
potential V to reproduce lattice-like data analogous to 
those in Fig. [U In the following we will call this method 
fit method. We will also discuss later the fit method for the 
present problem but first we want to propose a different 
strategy which does not require to assume a specific shape 
of the potential V, unlike the fit method. In the following 
we will call this new method direct method 20]. The 
idea of the direct method is to evaluate directly the irp 
amplitude using the expression 



T wp (E) 



1 



Gtt P {E) — G 7!p (E) 



(35) 



in the p max —> °o limit with Gk P (E) and G 7rp {E) from 
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FIG. 4. Synthesized lattice-like data considered in order to 
illustrate the inverse problem methods. Stable p case: (a). 
Unstable p case: (b). The error band represent the result of 
the fit in the fit method described in the text. 

Eqs. (|3U|) and for the energies of the points in Fig.Q] 
Recall that Eq. (|35|) is valid only for the energies solu- 
tion of Eq. ([M]) . Note that despite the fact that G vp (E) 
and G- Kp (E) are divergent in the limit p m ax — > oo, the 
difference G np (E) — G 7rp {E) which appears in Eq. 
is convergent. Therefore Eq. (|35|) is cutoff indepen- 
dent. This is definitely a non-trivial result an illustrates 
one of the strong points of the method to solve the in- 
verse problem. For practical numerical evaluations we 
have checked that considering an average in the interval 
Pmax ~ [1.5,2.5] GeV we get the same numerical result 
than considering p max — > oo. As shown in [20(, Eq. (155]) 
is a different and practical way of writing Liischer's equa- 
tion, taking into account the full relativistic two body 
propagators. 

In Figs. [5J and [5] we represent by square points with 
error bars the np amplitude (modulus squared, real part 
and imaginary part) for the stable (Fig. [5J) and unstable 
(Fig. [S]) p cases for the direct method. (The solid line 
and shadowed error bands represent the solution of the 




E [Me VI 

(c) 

FIG. 5. The up scattering amplitude solution of the inverse 
problem for the stable p case. Horizontal errors of ±10 MeV 
in the energy axis have to be also assumed. Square points 
with error bars: direct method. Solid line with error band: fit 
method 



fit method which will be explained later on). The central 
points of the solution of the direct method are obtained 
with the central values of the lattice-like set and the er- 
ror bars are evaluated by varying the energies within the 
10 MeV errors given in Fig. [U For the stable case there 
are no data points between 900 and 1100 MeV since, as 
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FIG. 6. Same as Fig. [5] but for unstable p meson 



can be seen in Fig. [31 one should go to very high values 
of the size of the box in order to get a good resolution in 
this energy region. It is worth noting that this method 
produces large errors for certain energies, particularly in 
the case of the unstable p. This is due to the fact that 
G- Kp {E) is usually very steep close to the energy values 
of the levels (see Fig. [2]) and, hence, small variations in 
E provides large variations in G„ p (E). This feature can 
be seen in the G^p vs. E plot in Fig. [5] where it is visible 



that the crossing between the V~ x and G vp lines usually 
occur close to poles of G^p and hence in an energy region 
where G wp changes rapidly. Anyway, a clear resonance 
shape corresponding to the ai(1260) is visible as a peak 
in |T| 2 . Note, however, that the shape is far from being a 
Breit-Wigner. The chiral unitary approach in which our 
lattice-like data set is generated, provides not only poles 
but the full scattering amplitude in the complex plane, in 
particular also in the real axis, and generates the possi- 
ble background besides the pole, i.e., provides the actual 
shape of the amplitude. Actually, for the ai(1260) there 
is a very strong background which distorts the shape from 
a Breit-Wigner. The reason is that the amplitude is zero 
at around E ~ M p since the potential V is zero around 
that energy, as can be easily seen just looking for zeroes 
of Eq. Q for the irp case. However, the pole contribu- 
tion itself has a large strength at this energy. This means 
that there must be a very strong background in order to 
cancel at that energy the tail of the Breit-Wigner shape 
coming from the pole in order to produce the zero. (See 
a more extended discussion on this zero in Ref. [44|). 

H0 [ J 



100 - 




E [MeV] 

(b) 

FIG. 7. Phase shifts obtained from the solution of the inverse 
problem for the stable (a) and unstable (b) cases. Horizontal 
errors of ±10 MeV in the energy axis have to be also assumed. 
Square points with error bars: direct method. Solid line with 
error band: fit method 
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From the scattering amplitudes we can get the phase 
shift, 5, that is well defined for stable particle an which 
in our normalization it is related to the amplitude via 



as 



„2i<5 



P 



T{E) = 



= 1 

AttE 
8irE 1 



■T 



P 



cot 5 



(36) 



where p = \J(E 2 — (m p + m v ) 2 ){E 2 — (m p — m n ) 2 ) / (2E) 
is the CM momentum. Note that Eq. (|3(j|) is only well 
defined for stable scattering particles. However, for the 
sake of comparison between the stable and unstable 
case we also use, by definition for the unstable case, the 
same equation using for the p meson mass the physical 
value, but being aware that in the unstable case this is 
not a well defined observable and it is just a theoretical 
mathematical exercise. The solutions of the inverse 
problem for the phase shifts using the direct method are 
shown by the squares and error bars in Fig. [7] for the 
stable, 7(a)| and unstable, |7(b)[ cases. (Again the solid 
line and shadowed error bands represent the solution of 
the fit method explained below.) 

It is interesting to note that the phase shifts that we 
have obtained are quite different from those of a Breit 
Wigner, where the phase shift would go from zero to 180 
degrees passing through 90 degrees at the pole of the res- 
onance. We see that the phase shifts stabilize around 
90-100 degrees at high energies around 1400 MeV. On 
the other hand, the shape of |T 2 | in Fig. Ois typical of a 
resonance, and bumps like that would be used to identify 
the resonance experimentally. To the light of this, and 
the comments above that the irp amplitude has a strong 
background at low energies, one might question the pro- 
cedure used in [3l[ where the shape of the amplitude is 
constructed from a calculated phase shift below thresh- 
old and another one above, assuming that one has a Breit 
Wigner amplitude. On the other hand, given the peculiar 
shape of the amplitude predicted here, the evaluation of 
phase shifts at physical energies in the energy range of 
our Fig. would be most welcome. In this sense, it is 
interesting to note that the approach of [31] produces S 
around 90 degrees in the region around 1430 MeV, where 
we find the phase shift stabilizing around that value (see 
Fig. [?}. 

Let us now apply the other method already introduced 
in Refs. [2(| to solve the inverse problem based on 
fitting the potential, fit method. 

The shape of the lowest order irp potential, V, based 
on Eq. ([3]) which comes essentially from chiral symmetry, 
can be written in the form [33j: 



V = a' + b's + -. 

s 



(37) 



In order to make the coefficients a, b and c adimensional 
and of natural order 1 we redefine the previous equation 



V 



-1 
4^ 



? a + b(s — m H ) -c 



(38) 



with mji = 1.2 GeV. Actually in Ref. 33] a, b and c were 
around 2, 1 and 3 respectively for the present channel, as 
can be obtained from Eq. Q. One should note that the 
expression of Eq. (|38p is the one that one has in the chiral 
unitary approach. In this sense, the fit that we perform 
provides as a best solution the results of the chiral unitary 
approach with Xmin = 0- O ne can gi ye one se lf some 
freedom to have other possible potentials. When this is 
done what one finds is that the best solution is essent ially 
the same but the uncertainties are somewhat larger [291 ]. 

Next, we assume that lattice data are provided by our 
synthetic lattice-like set, points with error bars in Fig. 2] 
Then we fit these points with the solutions coming from 
Eq. (f3~4"]) in order to get the best a, b and c parameters, 
which produce the minimum % 2 , Xmin- 

The G 7rp (E) function is of course dependent on the cut- 
off p m ax and, therefore, also are the parameters fitted. 
However, the amplitude T{E) obtained from Eq. |3"2"1) 
should be independent of this cutoff and therefore the 
inverse method does not require the knowledge or as- 
sumption of any particular cutoff, i.e., changing the cut- 
off value would produce different values for a, b and c 
but the amplitudes and observables derived from them 
would be the same. We have checked numerically that 
this is indeed the case within errors. This feature is not 
novel, it was already observed in [20] and is related to the 
behavior of amplitudes within the renormalization group 
method, as used for instance in Quantum Mechanics in 

In Fig. 2] we show with the solid line the result of the 
fit to the energy levels using for the fit the cutoff value 
Pmax = 1 GeV. The shadowed band represents the as- 
signed errors obtained by varying the fitted parameters 
such that x 2 < Xmin + 1- From Fig. till Fig. [7] we rep- 
resent by the solid line and error bands the results of the 
fit method for the irp amplitude and phase shift. Note 
that the fit method gives generally smaller errors than 
the direct method since it considers a global fit to all the 
points of the lattice-like set and does not suffer from the 
pathological error sensitivity discussed above for the di- 
rect method. However, it has the drawback of having to 
assume a shape of the potential, V, like in Eq. (|37J) . 

It is remarkable that in spite of the proliferation of 
eigenenergies in the case of the unstable p and the dif- 
ferent shapes of the two levels in Fig. 2] for the stable 
and unstable p, the amplitudes and the phase shifts in 
Figs. O [S] and [7| are rather similar. We observe then that 
the effects of considering the p unstable are rather mod- 
erate in the previous analysis. This is not the case if we 
start from data generated with unstable p but the inverse 
problem is analyzed with stable p, as will be explained at 
the end of this section. It is interesting to observe that 
the proliferation of levels in the case of the unstable p has 
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made the region of energies between threshold and 1400 
MeV accessible, while for the stable p the first accessible 
physical energies are around 1100 MeV. The price one 
pays for having these energies available in the unstable p 
case is a bigger uncertainty in the reconstruction of the 
phase shifts as seen in Fig. \7\ 

At this point let us discuss a different analysis which 
may be closer to what a lattice calculation implements 
(see, however, some technical observations at the end of 
this section). In a real lattice calculation the generated 
data should be closer to those in Fig. |4(b)[ unstable p 
case, since starting from quarks in the lattice and im- 
posing on them the boundary conditions would naturally 
lead to boundary conditions in all the pions, the spectator 
and those coming from the p decay. However, for simplic- 
ity, one may be tempted to do the subsequent analysis 
using the model for the stable p. In view of this, we are 
going to see what happens if one starts from a lattice-like 
set generated with the unstable model but the inverse 
problem is analyzed considering the stable case. Let us 
first focus on the direct method. Thus we evaluate now 
Eq. for the case of stable p-meson but for the en- 
ergies shown in Fig. |4(b)| The results for the amplitude 
and phase shift are shown in Fig. [8] The result is very 
far from reality. The real part of the amplitude blows-up 
for low energies and it has little to do with the results of 
Figs. [5] or O The phase shift is also different to Fig. 
even the sign. Conceptually there are no deep reasons 
why the results should be so different, but numerically 
the method is pathological for the following reasons. Let 
us focus on one energy where the extracted amplitude 
is very large, for instance the point of the lower energy 
ban d at L = 2.4 m; 1 , for which E = 839 MeV, (see 
Fig. |4(b)j ). Let us look now at Fig.[5J where we show the 
same plots as in Fig. [2 but for L — 2Am~ 1 , (Fig. [2 was 
evaluated for L = 2m; 1 ). The energy E = 839 MeV 
is indeed the first crossing point between V^ 1 (dashed 
line) and G^p (II ^ 0) (solid line). However, at this en- 
ergy the value of G 7rp (n = 0) (dashed dotted line) is very 
different, about a factor six smaller. In Eq. (|33|) . G np is 
similar both for stable and unstable case and has a small 
value at this energy, (about -0.02). Therefore the ampli- 
tude evaluated with G Tp (II = 0) is much larger than the 
one evaluated with G^p (II 7^ 0). On the other hand, the 
extra poles in G^p (II 7^ 0) coming from the eigenener- 
gies of the nn system as explained above, make the plot 
of G„ p (fl ^ 0) very different from G„ p (fl = 0). Thus, 
close to an extra pole this numerical analysis is also in- 
appropriate, since one would be misidentifying the poles 
responsible for the eigenvalues using the two procedures. 
These pathological behaviors of the G^p functions within 
the box are essentially due to the abruptness produced by 
the poles of the G function in the box and it makes very 
difficult to extract reliable information from this analy- 
sis. We have also implemented the fit method in the last 
analysis, with stable p but starting from data created 
with unstable p. In this case, no good fit is obtained 



since with this shape of the potential it is not possible 
to generate a curve decreasing with L as the lower level 
in Fig. |4(b)| does. Actually the best fit has ~ 200 
and produces the pole for the a\ in the real axis below 
threshold which makes also this analysis unreliable. The 
main conclusion of the analysis using the stable case but 
starting from data generated using unstable p meson is 
that the inverse method should be as close as possible 
to reality, which in this context means to consider the 
unstable p meson in the analysis of the data. Otherwise 
unreliable results could be obtained due to the abrupt 
shape of the G functions in the box. 

In practice, the above mentioned problem could not 
show up in some actual lattice calculations, like in [3l| . 
if large pion masses are used. Indeed, we have seen that 
with the pion mass = 266MeV, used in [3l|, the 
ground state level is much less affected than in the case 
shown here, and the extra level originating from II 7^ 0, 
which in Fig. [9] corresponds to the first excited state 
(eigenenergy around 1 GeV in the figure), is moved be- 
yond the level corresponding to the first excited pir state. 
In this case the results with the analysis with stable p 
would be similar to those with the unstable p. Further- 
more, one should also take into account that actual sim- 
ulations like the one of [31| . which do not incorporate 
three-pion correlators, would not see the decay of the p 
into 7T7T, in which case, the analysis with a stable p would 
be more appropriate to interprete such lattice results. 



VIII. SUMMARY 

We showed how to tackle the problem of the inter- 
action of two particles quantized in a finite box of size 
L when one of the them has a finite width. The idea 
is based on extending previously known techniques for 
the stable case but quantizing also inside this finite box 
the decay channel of the unstable particle. In this way, 
the continuous integration needed to evaluate the selfen- 
ergy is substituted by a discrete sum over the allowed 
levels in the box of the decay channel of the unstable 
scattering particle. We illustrate the method with the 
irp scattering which generates dynamically the ai(1260) 
resonance within the chiral unitary approach. The scat- 
tering energy levels inside the box for periodic boundary 
conditions, both for the stable and unstable case are com- 
pared. The results show significant differences between 
both cases. 

Then we explain how to solve the inverse problem 
of getting physical observables in the real world (infi- 
nite volume) from lattice data which are computed in 
a finite box. The idea is based on the improvement of 
the Liischer's approach developed in Ref. (20j, properly 
adapted to the present problem. We apply two methods 
to solve the inverse problem: the first one is a previously 
proposed way based on fitting the parameters of a given 
potential to get the lattice data energy levels, fit method. 
The second one, direct method, is to evaluate directly the 
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irp amplitude using Eq. (|35p . The advantage of the direct 
method, which is closer to the original Liischer approach, 
is that there is no need to provide a specific shape of the 
potential since there is no potential involved. However, 
the drawback is that the errors are larger than in the fit 
method and that the observables in the infinite volume 
can be evaluated only for those energies for which there 
are lattice data. 

With respect to the particular case studied here of the 
irp scattering and the a\ resonance, it is quite instructive 
to observe that the amplitudes obtained are rather pecu- 
liar and they do not resemble much the shape of a Breit 
Wigner. This feature will have to be considered in future 
lattice QCD studies. We showed that using the stable 



p, the first physical energies accessible were around 1100 
MeV. However, the proliferation of energy levels due to 
the quantization of the decay products of the p, making 
the study with the unstable p, made a wider range of 
energies available, although the induced phase shifts had 
larger errors using the direct Liischer method. In such a 
case, the fit method proposed here could provide a more 
efficient method to induce the phase shifts with much 
smaller errors. 

Furthermore, we have also discussed that if the start- 
ing generated lattice data takes into account the decay 
of the p meson into two pions, then the analysis of the 
inverse problem must be performed considering also the 
instability of the p meson. Otherwise the results could 
be numerically unreliable. 

The considerations done in the present work should 
spur the future lattice calculations to consider also the 
discretization of the decay channels when unstable par- 
ticles are involved in the scattering. In between, calcula- 
tions along the line of [3l| , which can be studied with the 
stable p analysis, produced at other volumes or within a 
moving frame, could bring extra information on phase 
shifts in the physical region which would eventually de- 
termine the peculiar shape of ixp amplitude in the a\ 
resonance region, with a strong diversion from a Breit 
Wigner shape. 
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